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Abstract — In this paper, we establish traveling wave 
solutions by using sine-cosine function method for Classical 
Boussinesq (CB), and the Mikhailov-Shabat (MS) equations 
which are nonlinear partial differential equations and also 
important soliton equations. It is shown that the sine-cosine 
method provides a powerful mathematical tool for solving many 
nonlinear partial differential equations in mathematical 
physics. 

Index Terms — Sine-Cosine function method, Traveling wave 
equation, Classical Boussinesq equation, Mikhailov-Shabat 
equations. 


I. INTRODUCTION 

The nonlinear evolution equations have a major role in 
various scientific and engineering fields, such as fluid 
mechanics, plasma physics, optical fibers, solid state physics, 
chemical physics and geochemistry. Nonlinear wave 
phenomena of dispersion, dissipation, diffusion, reaction and 
convection are very important in nonlinear wave equations. 
In recent years, quite a few methods for obtaining explicit 
traveling and solitary wave solutions of nonlinear evolution 
equations have been proposed. A variety of powerful 
methods, such as, tanh sech method {Malfliet W. [1], Khater 
et al. [2], and Wazwaz [3]}, extended tanh method {El-Wakil 
et al. [4], Fan [5], Wazwaz [6]}, hyperbolic function method 
{Xia and Zhang[7], and Yusufoglu and Bekir [8]}, Jacobi 
elliptic function expansion method {Inc and Ergut[9]}, 
F-expansion method {Zhang [10]}, and the First Integral 
method {Feng [11], Ding and Fi [12]} . The sine-cosine 
method {Mitchell [13], Parkes [14], and Khater [2]} has been 
used to solve different types of nonlinear systems of PDEs. 

The aim of this paper is to find the exact and periodic 
solutions by the sine-cosine method of the Mikhailov-Shabat 
(MS) equation and Classical Boussinesq (CB) equation, 
which are the important soliton equations. 

II. THE SINE-COSINE FUNCTION METHOD 

We introduce the wave variable £ = x - ct into the PDE 

V (u, U t U x u tt u KJ u tXi .. . ) = 0. (l) 

Where u(x,t) is traveling wave solution. This enables us to use 
the following changes: 
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One can immediately reduce the nonlinear PDE (1) into a 
nonlinear ODE 

Q (u f F 44 . .. ... ) 0 

The ordinary differential equation (3) is then integrated as 
long as all terms contain derivatives, where we neglect 
integration constants. 

The solutions of many nonlinear equations can be 
expressed in the form {Ali et al. [15], Wazwaz [16]-[17]} 

A shUfjuf j , 
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otherwise 
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Or, in the form 

pi co s ^GR) , 
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Where X, p and PtD are parameters that will be determined, 
p and c are the wave number and the wave speed respectively. 
We use 

u.('0 = A sin^(/i.*f) 

u n = X n silicon 

(u n ) £ = njipV cos(^)siti n|3-1 (^) (6) 

(u") ff = n 2 fj 2 p?A fl sin"* 1 Of) + 

n;i 2 i"/?(«/? - 1) sin^“ 2 Of) 

And 

■u(f) = /Lcos^Of) 
u n = 2" cos 00 

(u n ) £ = — n[i(3A n sinOO cos"' 3-1 (nf) 

(u n ') u - -n 2 ft 2 f? 2 r cos^Of) + 

nfi 2 A n p{?iP - 1) cos^“ 2 Of) 

And so on for other derivatives. 


We substitute (6) or (7) into the reduced equation obtained 
above in (3), balance the terms of the cosine functions when 
(7) is used, or balance the terms of the sine functions when (6) 
is used, and solving the resulting system of algebraic 
equations by using the computerized symbolic calculations 
We next collect all terms with same power in cos k ( pq) or 
sin^(p^) and set to zero their coefficients to get a system of 
algebraic equations’ among the unknowns p, P and X. We 
obtained all possible value of the parameters p, p, and X. 
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III. APPLICATIONS 

A. The Classical Boussinesq Equations 

The Classical Boussinesq (CB) equations 

u t + [(l + u)v]* = -\v xxx 

v t + vv x + = 0 (8) 

In order to obtain travelling wave solutions of equation (8), 
we make the transformations 
u(jc, t) = u(f); v(x, t) = v(£); if = x —wt 
And integrating once with respect to then (8) becomes; 


—wu + (l + ti)v +-v" = 0 

4 

And 


1 2 

u = wv - V 

2 


(9) 


( 10 ) 


Substituting Eq. (10) into Eq. (9) leads to the following ODE 

v" +4(1 — w 2 }v + 6 wv 2 — 2v° — 0 (11) 

Using (6) into (11) we get, 

—p 2 p 2 A sin? (fit;) +{rAp(p — 1) sin^ -2 (ft£) + 
4(1 — w 2 ')A sin^(prf) + 6wA 2 sin 2 ^(p^) — 
-2l 3 sin 3jS (juO = 0 

( 12 ) 

Equating the exponents and the coefficients of each pair of 
the sine functions we find the following algebraic system: 


(13) 


p((3-l) = 0 
p — 2 = 3p -* p = -1 

Substituting Eq. (13) into Eq. (12) to get: 

— irA sin _1 (jif) + 2 p 2 Asin~"(p^') + 4(1 — 
w 2 )2,sin -1 (/i£) + 6wX 2 sin _2 (j[/0 — 

2A 3 sin -3 (^if) — 0 

(14) 

Equating the exponents and the coefficients of each pair of 
the sine function, we obtain a system of algebraic equations: 


sin d (pC) : 2\irA — 27<r = 0 

—2/..SV &wA~ = 0 


jin 


sin 1 (|^0 : —p, 2 A + 4(1 — tv 2 )A = 0 

By solving the algebraic system (15), we get 


(15) 


3 — —1; p — +2y (l — w) ; A — 

+2-^/ (1 — w 2 ) 

(16) 

In view of (4) and (16) we obtain the periodic solutions 
v = +2-^ (1 — w 2 ) sin -1 (+2-^/(1 — w z ) A 

(17) 

Where f = i- wt 

To find the solutions u (x. t) , according to (10) 


1 

u = wv - V 

2 


2 


(IB) 


By means of the equations (4) and (16) and using equation 
(18), we have the following periodic solutions foru(:r, tj: 


u = +2q (1 — tv 2 ) w sin 1 (+ 
2(1 — w 2 ){sin -1 (pif)} 


(/■! _ 


2 


w 2 ) 0- 


(19) 


Where £ = x —wt. 

From the equations (17) and (19) we get the solutions of the 
CB equation by the sine-cosine method. 

B. The Mikhailov-Shabat (MS) Equation 

In this section we deal with the Mikhailov-Shabat (MS) 
equations 

Pt = Vxx + (p + ?)?* - + ^) 3 

1 

-?t = Vxx ~(p+ - - (p + ?) 3 

In order to solve MS system (20) we now introduce the 
transformation 


u(x, t) = p (x, t) + q (x r t) 
v{x, t) - q x (x, t) - p x (x, t) 
Then the MS system (20) becomes 


( 21 ) 


u t + v x — uu.. = 0 

v t + (m), - U 2 u x + u xxx - 0 (22) 

Substituting 

u(x f t) = = k(x — ct ) 

(23) 

(Where k and c are real constants) into (22) and integrating 
once with respect to £, we get 
1 ? 

v = cu |-u L 

2 (24) 

And 

1 


k^u" + uv -u 3 — cv — 0 

3 

Taking (24) into (25), we obtain the following ODE; 


(25) 


, 2 


u" + 


—— T^ u V — 
2k c k 2 6k 


u 3 = 0 


(26) 


Seeking solutions of the form (6) we get: 


— p 2 ft 2 A sinP (p£) + p 2 Ap(p— l)sin ff- -(^f) + 
^tA 2 sin^Of) - ^ Asm* (/if) + 
sin^fcif} = 0 


(27) 

Equating the exponents and the coefficients of each pair of 
the sine functions we find the following algebraic system: 

p(p -1) + o 

Jff - 2 = 3/? -» 0 = -1 (28) 
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Substituting Eq. (28) into Eq. (27) to get: 


~H 2 k sin + 2 jc Asm 3 (£tf) + 

sin“ 2 (^0 - ^Asin” 1 ^) + 

ib^ 3 sin_3 00 = 0 


v(x r £) — +v 12 c 2 esc [i J—b c 0] + 


6 c 2 esc"[+ 




(36) 


And 


Equating the exponents and the coefficients of each pair of 
the sine function, we obtain a system of algebraic equations: 


sin a (|i?) : 2|i a A + ^ 
sin“ 2 (|4> ^A 2 =0 

t«f 

sin -1 (|^) : —(J. : A — A = 0 


1 


A 3 = 0 


(29) 


By solving the algebraic system (29), we get, when — < 0, 


(3 — —1; jLt — +V12 c ; A — +Yl2 c 

"" 

In view of (4), (5), (23) and (30), for — < 0, we obtain the 
periodic solutions 


= +V12 ccsc[+ |— —r k(x — ct)] 


(31) 


Where, 0< +_ - — k(x — ct j <71 
And 


«(x,t) — +V12 csec[+ j — — k(x — ct)] 


(32) 


Where, 0< ± -77 k(x — ct') 

c z 

And for — > 0, we obtain the periodic solutions 


.2 


u{x,t) =+i\ 12 cesc h[+ — k{x — ct)] 


(33) 


And 


r 


u(x f t} = +Vl2 csech[+ — k(x — ct}] 


(34) 


To find the solutions v(x f t) , according to (24) 

1 . 

v = cu + — 

2 (35) 

By means of the equations (4), (5) and (23) and using 
equation (35), we have the following periodic solutions 
for v(x f t): 

c 2 

When — < 0.. we get 


v(x } i) = +Vl2 c 2 sec[±J—-^ k(x — ct}] + 


6 c 2 sec 2 [+J —— k(x — ct)] 


When — > 0, we get 


(37) 


v(x r t) = +ivl2 c 2 csch [+_J— — ct)] — 


6 c 2 csch 2 


±Jb k (x-ct) 


(38) 


And 


v{x f t) = +V12 c 2 sech[+J^r k(x — ct)] + 


6 c 2 sech 2 [+J— k(x — ct)] 


(39) 


IV. GRAPHICAL INTERPRETATION 

In this section, we will put forth the graphical 
representation of determined traveling wave solutions of the 
Classical Boussinesq (CB) equations. 
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Kb) 

Figure 1: (a) u(x,t) in equations (19) and (b) v(x,t) in 
equations (17) where vr = 0.25, a = 1 & taking +ve sign. 



2(a) 



2(b) 


Figure 2: (a) u(x,t) in equations (19) and (b) v(x,t) in 
equations (17) where vr = 0.25, a = 1 & taking -ve sign. 



3(b) 

Figure 3: (a) u(x,t) in equations (19) and (b) v(x,t) in 
equations (17) where u- = 0.3, p = 2 & taking +ve sign 



4(a) 
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4(b) 


Figure 4: (a) u(x,t) in equations (19) and (b) v(x,t) in 
equations (17) where vr = 0.3, p = 2 & taking -ve sign. 





6(b) 

Figure 6: (a) u(x,t) in equations (19) and (b) v(x,t) in 
equations (17) where w = 2 f u = 0.3 & taking -ve sign. 



In this section, we will put forth the graphical 
representation of determined traveling wave solutions of 
Mikhailov-Shabat (MS) equation 


*N_ 

□ 


1 


5(b) 

Figure 5: (a) u(x,t) in equations (19) and (b) v(x,t) in 
equations (17) where vr = 2, a = 0.3 & taking +ve sign. 
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7(b) 

Figure 7: (a) u(x,t) in equations (31) and (b) v(x,t) in 
equations (36) where k= 0.25,. c = 1. 


2 





9(b) 

Figure 9: (a) u(x,t) in equations (32) and (b) v(x,t) in 
equations (37) where k = 1.5, c = 3.5 


15. 


20 . 



Figure 8: (a) u(x,t) in equations (31) and (b) v(x,t) in 
equations (36) where k = 0.25.. c = —1. 
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10(b) 


Figure 10: (a) u(x,t) in equations (32) and (b) v(x,t) in 
equations (37) where k = 1.5, c = —3.5 



12(a) 




12(b) 

Figure 12: (a) u(x,t) in equations (33) and (b) v(x,t) in 
equations (38) where k = —2.5, c = —l.S 



Figure 11: (a) u(x,t) in equations (33) and (b) v(x,t) in 13( a ) 

equations (38) where k = 2.5, c = l.S 
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13(b) 


Figure 13:(a)u(x,t) in equations (34) and (b) v(x,t) in 
equations (39) where k = 4.2, c = 3.8 



14(b) 


Figure 14: (a) u(x,t) in equations (34) and (b) v(x,t) in 
equations (39) where k = 4.2, c = - 3.8 

V. CONCLUSION 

In this paper, we proposed the sine-cosine method to seek 
the periodic solutions of the NPDE. It is shown that this 
method is more powerful in giving more kinds of solutions. 
By making use of the method, we study the Mikhailov-Shabat 


(MS) equations’ and the classical system of Boussinesq 
equations’, new families of solutions written are found. The 
method is used to find a new exact and periodic solutions for 
the CB and MS equations. Thus, we can say that the proposed 
method can be extended to solve the problems of nonlinear 
partial differential equations which arising in the theory of 
solitons and other areas. These solution was similar to the 
solutions obtained in other paper. The study reveals the 
power of the method. 
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